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Abstract 



In this paper, weak symmetries of the Benjamin-Bona-Mahony (BBM) equation 
have been investigated. Indeed, this method has been performed by applying 
the non-classical symmetries of the BBM equation and the infinitesimal gen- 
erators of the classical symmetry algebra of the KdV equation as the starting 
C"| . constraints. Similarity reduced equations as well as some exact solutions of the 

BBM equation are obtained via this method. 
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>, 
CN , 1. Introduction 

^H , 

The Benjamin-Bona-Mahony equation 
p— ■ BBM : ut + u x + uu x - u xxt = 0, (1.1) 

o. 

used to model an approximation for surface water waves in a uniform channel 
|l| . If we note the KdV type equation 



u t + u x + uu x + u xxx = 0, (1.2) 

then we find out the likeness between these equations. Indeed, this similarity 
is not stochastic. Both of them used to model the waves appear in liquids, 
compressible fluids, cold plasma and enharmonic crystals which are of surface, 
hydro-magnetics, acoustic-gravity and acoustic types, respectively. The inter- 
esting point is that the main difference between equations (|l.ip and p. 21) occurs 
in the case of short waves (Find more information in |l|, |3j). 

The physical applications and mathematical properties of the BBM equa- 
tion (|l.lj) have been motivated many investigations such as obtaining the exact 
solutions via finite difference discrete process, global attractor and etc. 
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In this paper, we find the similarity reduced ODEs as well as resulted sim- 
ilarity solutions of this equation via weak symmetry implementation. Indeed, 
the organization of the present paper is as follows: Some historical information 
on the weak symmetry method are given in section 2. In section 3, we follow 
[10 ] in order to describe the theory of weak symmetries. Section 4 is devoted 
to performing this new class of symmetry methods using the invariant surface 
condition of the BBM equation (which is indeed the non-classical symmetry 
method) and infinitesimal generators of the classical symmetry algebra of the 
KdV type equation as the starting points in the weak symmetry method im- 
plementation. Finally, we have compared our results with those related papers 
using the classical symmetry method in order to clarify the advantages and 
disadvantages of the both strong and weak symmetry methods. 

2. Background 

Symmetry methods for differential equations, was originally developed by 
S. Lie [7j. These methods without any doubt are very useful and algorithmic for 
analyzing and solving linear and non-linear differential equations. Classification 
of differential equations as well as linearization of them are some other important 
applications of the symmetry transformation approach. First G.W. Bluman 
and J.D. Cole introduced the notion of the non-classical symmetry group of 
differential equations specially for the heat equation in 1969 (Find more infor- 
mation in (2|). For the non-classical method, we seek the invariance of both the 
original equation and its invariant surface condition, exactly this constraint (i.e 
invariance surface condition) causes the non-classical solutions which are more 
general than the classical ones. There are various implementations for perform- 
ing the non-classical symmetry method, for example, using the compatibility 
condition has been suggested by G. Cai and X. Ling [5J. 

First the weak symmetries have been introduced by P.J. Olver, and 
P. Rosenau in 1986 as a generalization of the non-classical symmetries with 
motivation of finding every solutions of the given system. In principle, not only 
the invariant solutions corresponding to arbitrary transformation groups can be 
found by the reduction method, but also every possible solution of the system 
can be found by using some transformation groups. In other words, there are 
no conditions that need to be placed on the transformation group in order to 
apply the basic reduction procedure (Find more information in (l0(). In the 
next section, we have an attempt to explain the notation and implementation 
of the weak symmetry method by considering the BBM equation as an example 
in order to prepare an appropriate setting. 

3. On the weak symmetry method 

Symmetry groups of a system of partial differential equations can be 
defined in two types (see [1C|). 



Definition. Let A be a system of partial differential equations. A strong sym- 
metry group of A is a group of transformations G on the space of independent 
and dependent variables which has the following two properties: 

a) The elements of G transform solutions of the system to other solutions of 
the system. 

b) The G— invariant solutions of the system are found from a reduced system 
of differential equations involving a fewer number of independent variables 
than the original system A. 

Definition. A weak symmetry group of the system A is a group of transfor- 
mations which satisfies the reduction property (b), but no longer transforms 
solutions to solutions. 

Indeed, there are several transformation groups which don't transform 
solutions of given equations again to solutions, but their differential invariants 
enable us to reduce them. In continuation we would illustrate the procedure 
of performing this method. For this purpose, first consider an arbitrary one- 
parameter transformation group, then substitute its related differential invari- 
ants and their derivatives into the original equation, finally, you will encounter 
with three different possible cases which in continuation have been illustrated for 
the (BBM) equation using an appropriate one-parameter transformation group. 

3.1. Reduced equation has no parametric variables 
Consider the one-parameter group 

{x,t, u) n> (x + A, t + A, u), 

So, we have the characteristic equation dx = dt = du/0. By substituting the 
resulted differential invariants i.e. r — x — t and w = u, into equation (|1.1[) . we 
have Wrrr + ww r — 0. As we see, this equation has no parameter variable. 

3.2. Reduced equation isn't incompatible and has parametric variables 
Consider the one-parameter group 

(x, t, u) n> (Xx, t, Am). 

So, the characteristic equation is dx/x = dt/0 = du/u. By substituting the 
resulted differential invariants r — t and w — u/x, in equation (jl.ll) . we have 
x(w 2 +w r )+w — 0, where w = is it's solution and this equation has x as the 
parametric variable. 



3.3. Reduced equation is incompatible and has parametric variables 
Consider the one-parameter group 

(x, t, u) H- (x + 2Xt + A 2 , t + A, u + 8Xt + AX 2 ). 

By substituting the resulted differential invariants i.e. r = x—t 2 and w = u— At 2 , 
in equation (|l.ip . we have ww r + w r + (8 — 2w rrr — 2w r )t + Aw r t 2 , where this 
equation has t as the parametric variable and it is incompatible. Indeed, from 
the coefficient of t 2 we have w r = and from the coefficient of t we have 
w rrr + w r = 4, this means that these equations are incompatible. 

4. Implementation of the weak symmetry method for the BBM equa- 
tion 

Since, the weak symmetry method is based on conjecture, so here, the several 
ideas of performing this method as well as some of its aspects are presented. 

4-1. Non-classical symmetries of the BBM equation 

There are several implementations to find the non-classical symmetries. 
Here, we follow the procedure presented by G. Cai et al. which they ob- 
tained the non-classical symmetries of the Burgers-Fisher equation based on 
the compatibility conditions [4|. 

Consider the following one-parameter group: 

x = x + e£(x,t, u) + 0(e 2 ), 

i = t + en(x,t,u) + 0(e 2 ), (4.3) 

u = u + £(p{x, t, u) + 0(s 2 ), 

Assume that the equation Ai(x,u^ n ') :— eq (|l.l[) is invariant under the trans- 
formation group (I4.3|) with the following invariant surface condition: 



A 2 (x,uW):=r)u t + Zu x -<p = (4.4) 

This means that ^ ( - 4 -'Ai|ai=o.A2=o = 0, where 

X = £(x, t, u)d x + rj(x, t, u)d t + ip(x, t, u)d u , 
is the infinitesimal generator of (|4.4[) . and 

x (4) =x + ( P x a Um + ... + <p mt a u 



J lHttt 1 



is the fourth prolongation of X, with the coefficients defined as tp J — DjQ + 
£,Uj x + nujt, where Q = (p — S,u x — rju t is the Lie characteristic and Dr = 
Si=o u Ji ®uj ^ ^ ne total derivative w.r.t. J (Find more information in [8j, |9() 

Without loss of generality in condition (|4.4[) . two cases n — and r\ = 1 
must be considered. 



Case I r\ = 1: In this case we have Ut = tp—£u x . Substituting this expression 
in (jl.ip we have D t ((p — S,u x ) = D t (u xx t — u x — uu x ), where D t is total derivative 
w.r.t. t. By substituting S,u xx in both sides of above, we find 

¥>' = Uxxtt - U xt - U t U x - UU xt + £,U XX 

= D xxt {u t ) -(u+ l)D x (u t ) + (£u x - (p)u x + £u xx , 

= D xxt ((p-£u x ) - (u + l)D x (cp - £u x ) + (£u x -ip)u x + £,u xx , (4.5) 

= p xxt - £u xxxt - (u + l)p x + {u + l)£u xx + (£u x - tp)u x + £u xx , 

By virtue of D x (u t ) = D x (u xxt —u x —uu x ), we have u xt = u xxxt -u xx —uu xx -u x . 
Finally, we find the following governing equation: 

^^-(un-iy -¥«*», (4.6) 

where tp 1 = D t (tp — £,u x ) + £u xt , tp x =D x (ip — £u x ) + £,u xx , and 

f xxt = D xxt {ip - £u x ) + £,u xxxt . 

By substituting the coefficient functions tp* , <p x , p xxt into invariance condition 
(|4.6[) . we are left with a polynomial equation involving the various derivatives 
of u(x, t) whose coefficients are certain derivatives of £ and (p. Since, £ and ip 
depend only on x, t, u we can equate the individual coefficients to zero, leading 
to these complete set of determining equations: £ x — £j = t; u — 0, ip = 0. So, 
we have £ = Ci, ip = 0. So, we find the infinitesimal generators of the non- 
classical symmetries using the above results as follows, when c\ = 1, we have 
01 = u x + Ut, and for c\ ^ the symmetries are 02 = u x , 0-3 = Ut- As a result 
we can state the following proposition: 

Proposition. The non-classical symmetries of the BBM equation in the case 
of t] = 1, spanned by 

01 = u x + u t , a 2 = u x , 0-3 = u t . (4.7) 



As a result of above proposition we have the following group-invariant solu- 
tions: 

1) For cti = u x + Ut, substituting it into u\ (u) we find u = F(x — t), where 
F must satisfy in: FF' - F'" = 

2) For (T2 = u x , substituting it into 02 (m) = we find u = F(t) for an 
arbitrary F, so from equation (|l.ip we obtain: u = 0. 

3) For 03 = Ut, substituting it into 0-3(11) = we find u = F(x), where from 
equation fTI]) F satisfies this equation: F' + FF' + F'" = 0. 



Case II rj = 0: In this case, without lose of generality we can let £ = 1, so we 
have: u x — (p. Using this we can deduce A(x,t,u) = <p x t — f — ucp. Subsisting 
this in the determining equation Atp u + (ft — A u (p — A x = 0, we obtain: 

ifixflfiu - 2ip(p u - Uiptfu + (f t = ipxtuf + uip u ip + 2tp 2 + ip xxt + ip x + wp x . (4.8) 

By assuming (p — tp(x, t) above equation changes into 

<p t - 2<p 2 - tp xxt -ip x - utp x = 0. 

So we have: ip = l/(c—2x). As a result, we deduce that u{x,t) = x/(c— 2t)+F(t) 
(where F is an arbitrary function) is a solution of (|1.1|) . 

If.,2. Using the classical symmetries of KdV type equation lil.2\) 

Since the appearance forms of equation (|1.1[) and (|1 -2|) are similar, we want 
to try our chance in order to obtain new similarity reduced ODEs for BBM 
equation through infinitesimal generators of the classical symmetries (CS) of 
KdV type equation as the starting constraint. For the classical symmetries of 
the KdV type equation using Lie classical symmetry we have the next theorem 
(Since, the proof is computational, to keep scope we don't present it here. Find 
more information in [alSl). 

Theorem. If we consider X = t;(x,t, u)d x + n(x,t,u)dt + i](x,t,u)d u as the 
infinitesimal generator of the classical symmetry group of the KdV type equation 
\1.2\l , then we have 

1 2 

n = c 1 t + c 2 , £ = -ci(x + 2i) +c 3 t + c 4 , <p = --c\u + c 3 , (4.9) 

where c\, ci, C3 and C4 are arbitrary constants. 
Hence the next corollary could be stated: 

Corollary. The classical symmetries of equation U.ty) i.e. KdV type equation, 
spanned by: 

X l = {x + 2t)d x +3td t -2ud u , X 2 = d t , X 3 =td x + d u , X± = d x . (4.10) 

So, we can consider any linear combinations of given vector fields in the 
above corollary as the starting constraint of the weak symmetry method. In 
continuation, we will illustrate the weak symmetry method using some linear 
combination of X\, X 2l X3 and X4 as the starting point. 

Example. Consider the one-parameter transformation group with the infinites- 
imal generator X 2 + X 3 — td x + d t + d u . The characteristic equation is dx/t = 
dt = du. So, we find the differential invariants as r = t 2 — 2x, w = u/t. 
By substituting these new variables in the original equation (jl.l[) we deduce 
(2w r — 2ww r — w rrr )t 2 — 2w r t + 4w rr + w, where t can be considered as the 
differential parameter. Note that solving the above ODE doesn't give new so- 
lution. 



Example. Consider the one-parameter transformation group with the infinites- 
imal generator X3 = td x + d u . The characteristic equation is dx/t = dt/0 = du. 
So, we can obtain the differential invariants as r = t, w = u — x/t. By substitut- 
ing these new variables in the original equation (11.11) we find: rw r -\-w — 1 = 0,, 
solving this reduced equation we obtain w = r/(r + c). So we can find u = 
(tx + x 2 + cx)/(t(x + c)) as the solution of equation (|1.1|) . 

4-3. Some other suggestions 

Some other ideas may be useful to reach other solutions of the BBM equation. 
For example, non-classical potential symmetry method or using classical and 
non-classical symmetries of other equations which have the similar forms as 
the BBM equation. Meanwhile, Physical knowledge of the model framework 
can be so effective in order to reach favorite solutions via weak symmetries. 
For example if you know your desired solution may be invariant under some 
scale of specific variables then the weak symmetry method can be started with 
an appropriate scaling transformation. Since the main goal of this paper was 
introducing weak symmetry method for BBM equation, we lay away performing 
of above approaches. 

5. More discussions 

Now, we want to compare our results with other related papers. Paper [6( 
is concentrated on the classical symmetries and optimal Lie system of the BBM 
equation. Comparing with [6|, we deduce that in this paper by applying the 
weak symmetry method we have obtained more similarity solutions and other 
useful suggestions are presented in order to reach more other solutions. 

Taking into account the sections 2 and 3 of [6f, the next theorem can be 
resulted (Find more information in ([8], Chapter 3). 

Theorem. If u = f(x,t) is solution of the BBM equation il.l]) . so are the 
functions 

u = f(x — e,t), u = f(x — as,t — s), u = e ( - u+1 ' )E f(x - as, e~ Et t), 

where e< 1 and a are arbitrary constants. 

Indeed, above theorem characterizes the invariant solutions of the BBM 
equation, for instance if u = c is a solution of equation (|1.1[) . then from this 
theorem we obtain u = ce 6 ^ u+1 ' as a solution of the BBM equation. For another 
example, if we consider the solution u — (tx + x 1 + cx)/(t(x + c)) of equation 
(|l.ll) . from this theorem we deduce that 

(t — e)(x — as) + (x — as) 2 + c(x — as) 
(t — e)(x — ae + c) 

(where e ^ 1 and a, c are arbitrary constants), is again a solution of BBM 
equation. By using such approach, we are enable to obtain more new solutions 
for the BBM equation. 



Conclusions 

In this paper, we have presented a comprehensive explanation of the weak 
symmetry method as the generalization of the classical Lie symmetry method. 
Indeed, we have performed the weak symmetry method for the BBM equation 
which has been fulfilled by applying the non-classical symmetries of the BBM 
equation and using the classical symmetries of the KdV type equation as the 
starting constraints. Also, the similarity reduced equations as well as some exact 
solutions of the BBM equation are obtained via this method. Finally, we have 
compared our results with papers using the classical symmetry method. Other 
suggestions for finding new exact solutions are also presented. 
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